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Abstract 

I The objective of the present work is to study a cosmological model for a spatially flat Universe whose 

04 ' constituents are a dark energy field and a matter field which includes baryons and dark matter. The constitu- 

^ ' ents are supposed to be in interaction and irreversible processes are taken into account through the inclusion 

^ , of a non-equilibrium pressure. The non-equilibrium pressure is considered to be proportional to the Hubble 

■ parameter within the framework of a first order thermodynamic theory. The dark energy and matter fields are 

coupled by their barotropic indexes, which are considered as functions of the ratio between their energy densities. 
The free parameters of the model are adjusted from the best fits of the Hubble parameter data. A comparison 
of the viscous model with the non-viscous one is performed. It is shown that the equality of the dark energy and 
matter density parameters and the decelerated-accelerated transition occur at earlier times when the irreversible 
processes are present. Furthermore, the density and deceleration parameters and the distance modulus have the 
correct behavior which is expected for a viable scenario of the present status of the Universe. 
; PACS:98.80.-k, 98.80.Jk 
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'■ 1 Introduction 

00 . ^ 

, The cosmic observations from the type la supernovae [T] suggest that the present period of the Universe is ex- 
' perimenting an accelerated expansion. Since matter contributes with attractive forces and positive pressure which 
decelerate the expansion, an exotic component - the so-called dark energy - with negative pressure must be postu- 
0^ ' lated to take into account the present accelerated expansion. Another dark component is also necessary to explain 
, the measurements of rotation curves of spiral galaxies [2]. This component, called dark matter, interacts only grav- 
itationally with ordinary matter. Dark energy can be modeled by a cosmological constant [3], however, it suffers 
the so-called fine-tuning and cosmic coincidence problems [3]. Hence, several models for the dark energy having 
dynamical properties were analyzed in the literature, among others we cite: scalar fields, tachyon fields, fermion 
fields, phantom fields, exotic equations of state and so on. 
;h ' It is expected that the dark components do not evolve separately. Indeed it is known that the problems stated 

5^ , above have a promising resolution if we take into account a dark energy-dark matter interaction. This interaction is 
supposed to be negligible at high red-shifts while it is preponderant at lower red-shifts. This may also alleviate the 
coincidence problem in the sense that it is possible to choose an appropriate form for the interaction term leading 
to a nearly constant ratio between the energy density of the matter field and the one of the dark energy at low 
red-shifts. Several cosmological models were proposed with interacting dark components, among others we quote 
the works given in the reference [S]. On the other hand, it is also understood that irreversible processes in the 
evolution of the Universe may also contribute significantly for the alleviation of the coincidence problem. In the 
Friedmann-Robertson- Walker metric this is effectively done through the introduction of a bulk viscosity associated 
to a non-equilibrium pressure, (see e.g. [5]). 

The aim of this work is to develop a cosmological model for a spatially flat Universe with interacting dark 
components where irreversible processes are considered. We follow [S] and couple the dark energy and matter fields 
by their barotropic indexes, which are considered as functions of the ratio between their energy densities. This 
is in contrast with most works in the literature, which directly consider an explicit form for the interaction term. 
Furthermore, we introduce a non-equilibrium pressure - within the framework of a first order thermodynamic theory 
- which is the responsible for the irreversible processes. 

The work is structured as follows. In section [2] the general features of the proposed model for the Universe - 
where dissipative effects are present and with interacting dark fluids - is discussed. The analysis of the cosmological 
constraints and cosmological solutions which follow from the model is the subject of section [31 Finally, in section |3] 
we present our conclusions. 
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2 Dissipative interacting dark fluids 



Let us consider a homogeneous, isotropic and spatially flat Universe described by the Friedmann-Robertson- Walker 
metric ds^ ~ dt^ — a{t)^{dx'^ + dy^ + dz^) where a{t) denotes the cosmic scale factor. Furthermore, let us consider 
a cosmological model where the Universe is modeled as a mixture of two constituents, namely, dark energy (de) 
and matter (m) which represents the baryons and the dark matter. In this model irreversible processes are taken 
into account by considering a non-equilibrium pressure and it is supposed to exist an energy transfer between the 
dark energy and the matter field. The Friedmann equation and the evolution equation for the total energy density 
P = Pm + Pdc read 

3i?^ = Pm + Pdc, (1) 

An + Pdc + 3iJ(pm + Pdc + Pm + Pdc + n?) = 0. (2) 

In the above equations H ^ a/a denotes the Hubble parameter, p = +pdc is the total equilibrium pressure and 
w stands for the non-equilibrium pressure, also known as the dynamic pressure. 

We follow [5] and decouple into two "effective conservation equations" , namely, 

Pm + 3-ff7^,Pm = 0, Pdc + 3i^7dcPdc = 0. (3) 

Above, it was introduced the effective barotropic indexes 7? (i = m, de) related by 

7m = 7. + ^^^ + -, (4) 

r Pm 

where r = Pm/pdc denotes the ratio between the energy densities and 7^ (i — m, de) represent constant barotropic 
indexes of the equations of state pi = (7^ — l)pi. This decoupling is motivated from the fact that we do not assume 
an explicit form for the interaction term between dark matter and dark energy. Rather we consider this interaction 
is intrinsically connected to their barotropic indexes. 

Again by following [9l we assume that the effective barotropic index of the dark energy is given by 

7dc = 7dc - F{r), (5) 

where F{r) is a function which depends only on the ratio of the energy densities r. The physical motivation for 
this choice is given by the interaction between the dark fluids. Indeed while 7^ and 7^^ give the influence of the 
interaction term in the field equations, F{r) accounts for the nature of this interaction. Since we are concerned 
with the coincidence problem, it is reasonable to suppose that F depends on the ratio r ~ Pm/ Pde- By taking into 
account the previous representation for 7^^ we can rewrite ([3]) as 

Pm + 3if7mPm = SH pdcF - 3Hzu, (6) 
pde + 3i?7dePdc = 'SHpdcF. (7) 

Within the framework of ordinary (first order or Eckart) thermodynamic theory the non-equilibrium pressure 
(see e.g. [6]) is proportional to the Hubble parameter H with proportionality factor identified with the coefficient 
of bulk viscosity 77, i.e., w = —3ijH. According to kinetic theory of relativistic gases (see e.g. [7]) the bulk viscosity 
is proportional to the temperature with an exponent that depends on the intermolecular forces, so that it is usual 
in cosmology to assume that r/ oc p™, where m is a positive constant. 

If we suppose that the coefficient of bulk viscosity is proportional to the square root of the total energy density 
~ V = Voy/P with 7/0 a- constant - the field equations are integrable and the expression for the effective barotropic 
indexes (|4]) become 

7fn = 7m + ^^^^^^ - ^ (1 + 7) 

We may infer from ([5]) and ([8]) that the effective barotropic indexes are functions only of the ratio between the 
energy densities. 

Now let us analyze the evolution equation for the ratio between the energy densities, which is given by 



-3HrT{r), (9) 



where T{r) denotes the expression 

J-(r) 



7m - 7dc + (^1 + J:) [F{r) - y/3m) 



(10) 



If we assume that a stationary state of the Universe is attained by a constant value of r = r^, this imphes that 
J-"(rs) = 0. Hence, the constant solutions will be stable if 



dF{r) 
dr 



>0, (11) 



so that we obtain from (|10p the inequality 

r, (1 + r,) (^^) - (Firs) - V^Vo) > 0, (12) 

\ / r—rs 

by taking into account that the barotropic indexes 7m and 7de are constants. From the inspection of (|12p we may 
infer that the simplest choice F = ^/Sijo fulfills the above inequality. This choice fulfills the stability condition. 
Moreover, the interaction term in the form SHXp^e, with A a constant and proportional to pdc, is consistent with 
the Le Chatelier-Braun principle of thermodynamics as it was shown by [5] • 

In order to determine the solutions of the field equations, we start by analyzing the evolution equation for the 
energy density of the dark energy (j31)2- According to the ansatz ^ and of the choice of F, the effective barotropic 
index 7^^ = 7^0 — V^rjQ is a constant. Hence, we may integrate equation ^2 and obtain 

Pdc = Pdc [— ) > (13) 

where the index stands for the present values of the variables. 

From the differentiation of the Friedmann equation ^ with respect to time it follows 



The integration of the above equation leads to 



H + |(7,n - V3vo)H' ^(7.„ - 7dc)pSc [^) " = 0. (14) 



The constant of integration C is found by considering the current values of the cosmic scale factor ao and of the 
Hubble constant Hq, yielding 

C^H^- ^. (16) 
Hence, (|15p can be rewritten in terms of the red-shift z — (ag/a — 1) as 



El 

Hi 



(1 + + f}0^ (1 + z)37lo ^ (17) 



where 11^ = Pil{p-ai + Pdc) denote the density parameters. 

From the knowledge of — (pm + Pdc) /3, we can obtain the density parameters of the matter and dark energy 
in terms of the red-shift, namely, 

f7SJ,(H-z)3('rm-N/3r,„) +f^0^(l_^^)37=„' 

= ^JsilLLll . (19) 

170,(1 + z)3('^»-^^')») + 170^(1 + z)3^dc 

The determination of the the ratio between the energy densities as function of the red-shift follows from r{z) = 
^m{z) /^dc{z) ■ Furthermore, since the non-equilibrium pressure is given by nj = —3\/3riQH^ , it can be also expressed 
as a function of the red-shift thanks to P7|) . 

Another parameter which is important in cosmology is the deceleration parameter g = 1/2 -I- 3we/2, which is 
given in terms of the effective parameter We = (pm + Pde + 'cc)/{pni + Pdc)- From the barotropic equations of state 
and from the representation of the non-equilibrium pressure the effective parameter becomes 

We = (7m - mmiz) + (7^0 + ^VO - i)^dciz) - VS^lo- (20) 

By inspecting the expressions (|17p - (j20p we may infer that there exist three free parameters in the proposed 
model, which are the coefficients 7,„, 770 and 7^^,. In the next section an analysis to set cosmological constraints on 
the free parameters is performed and the cosmological solutions are analyzed. 



Table 1: Hubble parameter H{z) from [llj . 
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3 Cosmological constraints and cosmological solutions 

The coefficients 7m, '7o ^md 7^^, can be found from the observational cosmological constraints which are based on the 
data of the Hubble parameter li(z) given in Table [1]- taken from [11] - together with the values iJo = 72 km/(s 
Mpc), = 0.30 and = 0.70 [T^]. The set of values given in Table [T] was used in the work [5] and the adopted 
methodology is explained in the appendix. 

For the viscous case, we have considered a dust- like matter field (7™ = 1) and adjusted the parameters 7^^ and 
770. In Figure [T] it is plotted the probability ellipsis in the plane 7^^ versus 770 and the best fit value is indicated by 
a dot, which corresponds to 7^^ = 0.125445 and t]q 0.0140124 with = 9.104007. 

In order to interpret the results for the viscous case, we compare it with the non- viscous one, which refers also to 
a non-interacting model. In this case the free parameters are 7m and 7^^,, and in Figure [2] we show the probability 
ellipsis in the plane 7^^, versus 7m. The best fit value is indicated by a dot, which corresponds to 7!^, = 0.0259052 
versus 7m = 1.0051 with ^ = 9.1407510. 

In Figures [1] and [2] the points inside the inner ellipses or between them stand for the true values of parameters 
with 68.3% and 95.4% which correspond to \u and 2(7 confidence regions, respectively. 

From the knowledge of free parameters of the model it is possible to perform an analysis of the cosmological 
solutions. We start with the investigation of the density parameters which are plotted as functions of the red-shift in 
Figure[31 the solid lines corresponding to the viscous case, whereas the dashed lines to the non- viscous case. One can 
infer from this figure that the energy transfer from the dark energy to the matter field is more pronounced for the 
non-viscous case, since for this case the growth of the density parameter of the matter field and the corresponding 
decay of the dark energy with the red-shift are more pronounced than those for the viscous case. This behavior is 
expected when we analyze the evolution equations for the energy densities ([6]) and ([7]) for the viscous case and can 
also be verified from the analysis of Figure |4] which represents the evolution of the ratio of the two energy densities 
f = p-m/ PAb with the red-shift. This last figure also shows that in the future - i.e., for negative values of the red-shift 
- there is no difference between the two cases, since both tend to small values, indicating a predominance of the 
dark energy in the future. 

In Figure [S] the deceleration parameter is plotted for the two cases. The present values of the deceleration 
parameter g(0) and the value for the red-shift Zt where the transition from a decelerated to an accelerated regime 
occur are : (i) q{Q) ~ —0.43 and zt ~ 0.74 for the viscous case and (ii) g(0) ~ —0.52 and zt ~ 0.67 for the non- viscous 
case. These values are of the same order of magnitude of the values given in the literature: g(0) — —0.46 ± 0.13 
(see [13]) and zt = 0.74 ±0.18 (see [H]). When the viscous case is compared with the non-viscous one we may infer 
from this figure that: (i) the former has a smaller deceleration in the past than the latter; (ii) the transition from 
a decelerated to an accelerated regime occurs earlier for the former and (iii) in the present and in the future the 
former has a smaller acceleration than the latter. 

The effective index We as function of the red-shift is shown in Figure [S] for the viscous and non- viscous cases. 
We may conclude that in the future the mixture of matter and dark energy behaves like a quintessence and the 
non- viscous case approximates to a cosmological constant with We ~ —1. For large values of the red-shift tends 
to zero for the non-viscous case and to a small negative value for the viscous one. Here we call attention that for 
large values of the red-shift it is necessary to include a radiation field which will imply in a positive effective index. 

In Figure [7] it is represented the distance modulus /io, which is the difference between the apparent magnitude 
m and the absolute magnitude M of a source. Its expression is given by 

Aio-m-M = 51og|(l + z)^ _^|+25, (21) 
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Figure 1 



Confidence regions for tire best fit values for the viscous case. 



where the quantity within the braces represents the luminosity distance in Mpc. The circles in this figure are 
observational data for super- novae of type la taken from the work |15j . This reference contains 4 different data 
sets related to various light-curve fitters. For practical purposes, we adopted the SALT data set {Ry = 3.1). It 
is possible to conclude that there is a good fitting of the curve with the observational data. Moreover, it can be 
seen from the small frame in this figure that there is no sensible difference between the curves for the viscous and 
non-viscous cases. 

4 Conclusions 

In this work we studied a cosmological model with interacting dark fiuids in a dissipativc Universe where the non- 
equilibrium pressure is the responsible for the irreversible processes. The non-equilibrium pressure was supposed 
to be proportional to the Hubble parameter within the framework of a first order thermodynamic theory. The 
coupling between matter and dark energy was made through their barotropic indexes, which were considered as 
functions of the ratio between their energy densities. The function of the ratio between the energy densities - which 
is the responsible for the energy transfer between matter and dark energy - follows from the stability analysis of the 
differential equation for the density ratio. A procedure was performed to set observational constraints on the free 
parameters of the model by using the observational data of the Hubble parameter. It was shown that the energy 
transfer from the dark energy to the matter field is more efficient for the non-viscous case. Furthermore, for both 
the viscous and non- viscous cases we obtained that the dark energy density predominates in the future, the mixture 
behaves like a quintessence in the future and the values of the deceleration parameter are of the same order as those 
given in the literature. It was shown also that the behavior of the distance modulus fio ^ which is related with the 
luminosity distance - has a good fit with the observational values. 
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Appendix: Bayesian Inference 



In a statistical sense a physical model may be thought as described by a set of parameters. The determination 
of these parameters may be carried out in many ways; the most commonly used framework to accomplish this is 
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Figure 2: Confidence regions for the best fit values for the non- viscous case. 





Figure 4: Ratio between dark matter and dark energy as function of the red-shift z. Solid lines - viscous; dashed 
lines - non- viscous. 




Bayesian inference, a well-known method of statistical inference which employs evidence to estimate parameters of 
a model. The main purpose of this section is just to give a brief introduction to the subject. 

For a given model and data set, Bayesian inference employs a probability distribution called posterior probability 
to summarize all uncertainty. This probability distribution is proportional to a prior probability distribution (or 
simply the prior) and a likelihood function. The later, denoted by V(D\0), is usually defined as the unnormalized 
probability density of measuring the data D ~ {Di, D2, for a given model A4 in terms of its parameters 

9 = {01, 62, 0^}. For our purposes it suffices to assume that the measured values arc normally distributed around 
their true value, so that 

P(D|0)(xcxp[-x'W/2] . (22) 
The posterior 'P{d\'D) is determined by Bayes' theorem 

where V{9) denotes the prior probability distribution. The prior carries all previous knowledge about the parameters 
before the measurements have been performed. 

Parameter estimation is performed in Bayesian inference by maximizing the posterior V{0\D). This is in contrast 
with the frequentist approach, in which the likelihood 'P{T)\9) is maximized. Nevertheless, whenever the so-called 
uninformative priors are considered, both frameworks lead to the same conclusions. If the measured data are 
independent from each other as well as Gaussian distributed around their true value, D(0), then maximizing the 
likelihood V(D\9) is equivalent to minimize the chi-square function 

X^i9) = (D"'" - n{9))C-\T>°'" - T>{9)f, (24) 

where C is the covariance matrix given by the experimental errors. For uncorrelated data dj ~ Sijaf and 

where ai denotes the experimental errors. 

In Bayesian inference, the confidence intervals are drawn around the maximal likelihood point, giving the best 
fit parameters. It is conventionally used Icr and 2a confidence regions with 68,3% and 95,4% of probability, 
respectively, for the true value of parameters. These regions are mathematically defined by the inequalities 

x\d)-xH0,f)<2.3, (26) 

for 1(7 range and 

xHO)-x'{Obf)< 6.17, (27) 
for 2a range, where 9bf denotes the best fit value of parameters. 
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